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Path Antenna Gain in an Exponential Atmosphere 
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The problem of determining path antenna gain is treated here in greater detail than 
previously [1, 2, 3]. 1 The method used here takes into account for the first time the exponen- 
tial decrease of the gradient of refractive index with height, and a scattering cross section 
inversely proportional to the fifth power of the scattering angle. Results are given for all 
combinations of beamwidths and path geometry, assuming that symmetrical beams are 
used on both ends of the path and that atmospheric turbulence is isotropic. The result 
appears as a function of both of the beamwidths, in addition to other parameters, and thus 
the loss in gain cannot be determined independently for the transmitting and receiving 
antennas. The values of the loss in gain are generally lower than the previous estimates 
for which a comparison is possible. 

1. Introduction 

Theories of the forward scatter of radio waves in the troposphere or stratosphere usually 
differ only in their assumptions about the radio wave scattering cross section and pertinent 
meteorological parameters. With any theory, the total effective gain of highly directive 
antennas is less over a scatter circuit than it would be in free space. The "path antenna gain" 
is less than or equal to the sum of the free space gains G t and G r of the transmitting and receiv- 
ing antennas. Expressed in decibels above the gain of an isotropic radiation [1], 

G p <G t +G r . (1-1) 

Booker and de Bettencourt [2], Norton, Rice, and Vogler [1], and more recently Staras [3] 
have estimated the loss in antenna gain G L , where G L in decibels is given by 

(*l — Gt -h'.G r — G p , (1-2) 

for an atmosphere in which the variance of the gradient of the radio refractive index is either 
constant or varies parabolically with height. This paper estimates G L with a model atmos- 
phere in which the refractive index decreases exponentially with height. G L is evaluated using 
antenna radiation patterns which are the same horizontally and vertically. To the authors' 
knowledge, G L has not been evaluated using antenna patterns which are not symmetric about 
the beam axis. Figures 6 and 7 show graphs of G L for various values of the parameters, and 
section 4 summarizes important definitions and gives a simplified method for predicting G L . 
The reader who desires only an estimate of G L may skip to section 4. 

2. Antenna Gain in Free Space 

The gain function J(d, 4>) is the ratio of the power radiated in the direction (6, <j>) per unit 
solid angle to the average power radiated per unit solid angle. This expresses the increase in 
power radiated in a given direction by the antenna over that from an isotropic radiator emitting 
the same total power. It is convenient to express the gain pattern of an antenna and its free 
space gain in terms of the electric field distribution F(£, rj) over the antenna aperture. Figure 1 
shows a coordinate system with an antenna aperture in the x, y plane and centered at the origin. 



Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Geometry for determining the gain of an antenna. 



The point P(x, y, z) is in free spaee in the Fraunhofer or far zone. Using this notation, the 
power radiated per unit solid angle in a given direction by an antenna having constant phase 
distribution over the aperture is approximately [4] 



P(e,ct>y- 



/€//* 



2X 2 



J> 



77) exp [jksmd(£ cos 0+77 sin <£)]d£r/r7 



(2-1) 



where V/V € i s the impedance of the propagating medium, X is the free space radio wavelength, 
k=2ir/\, and d^d-r] is a surface element of the aperture. The average power radiated per 
steradian in terms of the total power P & across the aperture is 



PJ4*=dJ£fjF(i,r,)\»d&,, 



and the ratio of (2-1) to (2-2) is the gain function. 

M4,),=4*P(6,<l>)/P a . 



(2-2) 



(2-3) 



If the aperture illumination F(£,ti) is symmetrical about the origin, the maximum value 
/max(0,<£) of (2-3) corresponds to the direction 0=0 along the antenna beam axis, the 2-axis in 
figure 1, 

/mai(W) =/(<>,*) =^l,/(X 2 /47r) . (2-4) 

The factor A e , called the "effective absorbing area" of an antenna, is given in terms of the 
distribution function as 



^=ij/tt,ii)^ii|y£i^tt,ii)iw««ii. 



(2-5) 



For an isotropic radiator, A e =\ 2 /4:7r, as follows from (2-1) and (2-2). For any antenna with 
a uniformly illuminated aperture (F(%,ri) constant), the effective area equals the actual area, 
as may readily be seen from (2-5). For the more usual 10 db tapered illumination of the 
aperture of a parabolic dish antenna, the effective area A e is about 56 percent of the actual area. 
The free-space gain of an antenna in decibels above the gain of an isotropic radiator is 



(7=10 log 1( / max (0,*) = lO log 10 UV(X 2 /47r)]db. 
Define the normalized gain function, g m {B,<ji) : 



(2-6) 



(2-7) 
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Figure 2 shows the normalized gain function g m (6,<j>) versus 6 for a uniformly illuminated 
circular apert ure of 100 wavelengths diameter. Figure 2 would be the same for'any value of <£. 
Figure 2 also shows a "square" antenna power pattern, where g m (Bfi) maintains the value 
unitv for O<0<5 and is zero wiien exceeds 5. 
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Figure 2. Computed polar pattern {without side lobes) for a parabolic dish (circular aperture uniformly 

illuminated) with D/X= 100 

The "square" power pattern will be used throughout the calculations in this paper. 

For practical purposes the gain is often specified in terms of the half-power beamwidth, 25. 
For an antenna with a symmetrical distribution of field intensity over the aperture one obtains 
[4] 



25 = 



(2-8) 



whe re /z f is the ith moment of the function F(x) which gives the aperture field in terms of 
x=^e+y 2 - Combining (2-1), (2-7), and (2-8). 



or 



/(0^)-(mo/m 2 )/(25) 2 
jf(0 J *)=Cl/(2*)» J 



(2-9a) 
(2-9b) 



where d is a constant of proportionality depending on the characteristics of the antenna used. 
As examples, with the beamwidth 25 expressed in radians, Ci = 7 .2 for a half -wave dipole, and 
(7j =8.2 for a parabolic dish with the edge illumination 10 db down from the center illumination. 
This constant C\ can be determined for some types of antennas using tables in the book, Refer- 
ence Data for Radio Engineers [5]. 

The half-power beamwidth in figure 2 is 25, where 5 is the value of 6 corresponding to g m ($) = 
0.5 For the square pattern, assuming that 5 is small, the power is radiated through 7r5 2 stera- 
dians if the aperture A is circular, and through (25) 2 steradians if A is square. The power gain 
/(O,0) is then 4t/tt5 2 for the circular aperture or 47r/(25) 2 for the square aperture, that is, ^ = 16 
or Ci=4tt. 
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G Ly as determined by calculations in this report, depends directly on the semibeam widths 8 t 
and 5 r , and to a lesser extent on the shapes of the associated power patterns. Subscripts t and r 
in this paper refer to transmitting and receiving antennas, respectively. 

3. Scatter Integral 

If we denote by K the received power relative to free space of a given scatter propagation 
circuit with arbitrary antennas, and by K the received power relative to free space of the same 
circuit using isotropic antennas, then the loss in gain in decibels may be defined as 



o L 



-10 log 



K 
10 Ko 



and the total path antenna gain is 



Gp=G t +G-G L 



(3-1) 



(1-2) 



The loss in gain G L depends on the plane wave power patterns of the two antennas and upon 
the path parameters in such a way that the effect of one antenna cannot be evaluated inde- 
pendently of the power pattern of the other. 

In this report, it will be assumed that the antennas are sufficiently high so that ground 
reflection effectively doubles the power incident on the scattering volume over what it would 
be in free space and again doubles the power at the receiver. In general, K is given by the 
volume integral 



K-- 



where d is the straight line distance between transmitter and receiver, R and B are distances 
from transmitter and receiver, respectively, to a scattering element dv, and g t and g r are 2g mt 
and 2g mr , that is, double the power patterns defined by (2-7) . The scattering cross section a 
gives the fractional part of the incident power scattered in a given direction (at an angle, a-\-fi 
from the forward direction) per unit volume per unit solid angle. Refer to figure 3, which shows 




Figure 3. Geometry for the scattering integral. 
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some of the geometry appropriate for the evaluation of (3-2). In this figure and throughout 
(lie remainder of this paper, 6 is the angular distance, or the minimum scattering angle. A 
del ailed discussion of the derivation of the form for a used here may be found in a recent paper 
[6]. The following result is obtained: 



3X<(Ati) 2 > 

0= T-T^ 



[2sin(^)]- 5 , (3-3) 



where a+/3 is the angle between the directions R and Ro as shown in figure 3, X is the wave- 
length, <(An) 2 > is the variance of refractive index fluctuations in the scattering medium, 
and /o is a "scale of turbulence" appropriate for the scattering medium, assuming homogeneous 
and isotropic turbulence. Following Booker and Gordon [7], Abbott [8] derived in a sub- 
stantially different fashion the following expression for <(An) 2 >// 2 : 

<(An) 2 >/Z 2 = S / exp (-2 7 h)=S exp (-2 7 z) (3-4) 

where S f is the value of this parameter at sea level, h^z—h /2, S^S' exp (7A0), and 1/7 is the 
scale height of the atmosphere, here assumed constant. 

Using the coordinate system shown in figure 3, (3-2) becomes 



- J Xi J„ J-. [ 2 an (2±£)J(fifl»» (3 " 5) 

Using the approximation 

1 (dt-xyjda+x)* 

[2sin(^)]W" *0*+*>" ' (3 " 6) 

which is justified in appendix 1, and assuming that the effects of back scatter are negligible so 
that we may limit our x integration to the region between the transmitter and receiver, and 
limiting z to the region above the horizon rays, we obtain 



where 



and 



•>\ f r&i r™ r m ro r m / ,o ° "\ 

*-m.i. *L*L* M v-** , )J , L* M r (3 " 7) 

Zim^Po(d 2 +x) (3-9) 

22mSa (rfl — x) (3-10) 



s=% (3-11) 



We can express K in terms of four parameters, g mt , g, 

"a 

and 

v =t!0, (3-12) 

where 

= a„+/? o , (3-13) 

Thus 

K=K(s, v , g„ g r ). (3-14) 
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For isotropic antennas g mt =g mr =l, and thus g t =g r =2, including the effects of ground reflection 
as explained previously. K then becomes 

(3-15) is evaluated exactly in appendix 2 with the result given by (7-14). K (3-7) was eval- 
uated by imposing geometric limits on the integrals and using Ta}dor series expansions. The 
limits vary as su<^l, sn=l 3 or s/z>l where 



M= 



5 t 



(3-16) 



Figure 4 is a composite drawing and figure 5 is a drawing showing the xw and xh planes for 
,S'/i<1. This leads to five integrals corresponding to the regions shown in figure 5 since it 
was found necessary to evaluate the integrals corresponding to the regions I 42 and I 43 as one 
integral, and then subtract the integral for I 43 . The integral K was also evaluated numerical^ 
with <^(An) 2 ^>/ll proportional to z~ 2 and sr 4 , using the antenna power pattern [2J x (u)Ju] 2 of 
figure 2 (without side lobes), and the value obtained for G L was found to be always within one 
and a half decibels of the value obtained using a "square" pattern with the same half -power 
beamwidth and height dependence. 




Figure 4. Scattering volume with "square" antenna 
patterns. 




Horizontal Plane 



Figure 5. Shape of scattering volume in vertical and horizontal 
planes, sju>1. 



Effects of anisotropic turbulence have not been evaluated in this paper, but it can be shown 
that anisotropic turbulence arising from scale length changes has no effect on K Q [isotropic 
antennas] except to add a multiplicative factor of r } where 



l v _ scale of turbulence in vertical dimension 
l w scale of turbulence in horizontal dimension 



(3-17) 



However, it has been shown [6, 9] that anisotropic turbulence produces the following effect in 
the integrand, (3-8) of the scatter integral (3-7) 



fA{v)=g nt g, 



(vflJ+zH, 2 ) 6 ' 2 
g ml g m r (dr-xy^+xySe- 2 " 



/, 



(("M" 



(3-18) 
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where /^ (v) is/ (v) with the effects of anisotropy included. Using the square pattern as was 
done in (his paper, with geometric limits acting as a cutoff for the antenna beamwidths, it 
will be shown that, at least qualitatively, anisotropy will have some effect when using high 
gain antennas. As can be seen using (3-18), only the w integration is affected. A typical w 
integral using the geometric limits for anisotropic turbulence would be of the form 

p tf2+ *> *° (3 _ 19) 

w 
Let— =v. Then (3-19) becomes 

r 



/.' WfTW <**» 



This integrand is the same as in the case of isotropic turbulence, but the range of integration 
is larger, since r<l. Thus, the values of L obtained here would be lower if the effects of 
anisotropy were included. This effect is the same as if antennas with larger beam widths in 
the horizontal dimension than in the vertical dimension were used in an isotropic atmosphere. 

4. Results 

The integrals in the preceding section were calculated for numerous combinations of the 
parameters. Typical curves which resulted from these calculations are plotted in figure 6. 
The curves in figure 6 corresponding to path symmetry (s=l), and identical antennas (/x=l), 
are shown again in figure 7 using a different coordinate system which facilitates interpolation. 
Because of the number of parameters involved it was found impracl Leal to show all of the curves 
of the type shown in figure 6, and instead, an empirical formula was deiived, using these curves, 
which gives good accuracy using only two graphs and a few T simple calculations. The following 
discussion defines the necessary parameters and shows how to use the empirical curves. 

The basic parameters which depend on the path geometry are shown in figure 3 and are 
as follows : 

«o, ft, d, (4-1) 

from which may be calculated 

•=», (4-2a) 

d = a + Po (4-2b) 



and 

Combining these with the rate of decay of the variance of the refractive index with height, y, 
one obtains 

ri=ydd (4-3a) 

and 

v = 2yho= TT+^^ (4_3b) 

For the purpose of illustration, 2 y=0.23/mile will be used here. 



- For additional values of y, see reference [10]. 
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Figure 6. Loss in path antenna gain. 



50 100 




05 12 5 10 20 50 100 

Ratio 9/Zh of angular distance 9 to 
half-power antenna beamwidth, 28 

Figure 7. Loss in antenna gain, Gl assuming equal free space 
gains G t and G r at the terminals of a symmetrical path 5 t = 

5r = 5, s=l. 



The parameters which depend on the antenna characteristics are the half power semibeam- 
widths, 8 t and 8 r . These are used to calculate 



H=d T /8 t . 



(4-4) 



At the end of this section, formulas are given for calculating 8 from the free space gain of the 
antenna or from the frequency and area of the antenna aperture. 

In the empirical formula, the two parameters s and /x, are combined in the single parameter 
six, and two separate cases arise. 

Case I, sjjl>1. 

Calculate 



n=otol8 t . 



Read / (V) from figure 8 and calculate 



A=(n+Omv)/f(v). 



(4-5) 



(4-6) 



Read G L from figure 9 for the appropriate value of sp. Linear interpolation is possible using 
the formula 



G L =G L '+ y M j {0 L "-G L '), 



(4-7) 
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Figure 8. The contraction factor iiy) . 
Supplementary Values. 

v /O) » /(„) 

1 1.22 5 1.88 

2 1.41 15 3.10 
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Figure 9. Loss in path antenna gain versus d and the product 
sm for v = 0. 

where G L ' and G L " correspond to the graphically obtained values for (sp)' and {sp)" loss than 
sp and greater than sp, respectively. 

Case II, s/z<C 1 . 

Calculate 

Read/(V) from figure 8 and calculate 



n=p /5 r . 



(4-8) 



= (t*+0.03iO//(iO. 



(4-9) 



Read 6? L from figure 9 on the curve corresponding to the value of 1/sp. Use formula (4-7) for 
interpolating replacing (sp) by (1/sp), (sp)'/by (l/sp)' and (s/x)" by (l/sp)". 

If the gain ^ in decibels is known, the half-power semibeanrwidth may be calculated using 
the following formulas which may be derived from the formulas in section 2, and are given here 
for a parabolic dish of diameter D feet at a frequency j mc megacycles per second. G is given by 



Also, 



with h in radians, and thus, 



(9=10 log 10 [0.567r 2 (#/A) 2 ] 
= 20 log 10 Z^ + 20 log 10 / mc -52.431 db. 

<7=10log 10 [8.2/(25) 2 ] 

<5=1432X10- Gr / 20 =609X/£> mr. 
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(4-10a) 



(4-1 0b) 



(4-10c) 



5. Comparison with Previous Estimates 

In figure 10, our estimates of G L are plotted together with the estimates given by the previous 
authors [1, 2, 3]. The values from Staras' [3] paper which are shown in figure 10 correspond 
to the assumption of isotropic turbulence. For this comparison, it was necessar} 7 to use only 
the curves corresponding to the assumptions of path s^nnmetry and identical antennas on each 
end of the path since some of the previous estimates are available only for these cases. In 
addition, our assumption that the gradient of refractivity of the atmosphere decreases expo- 
nentially with height leads to the parameter v =2yh which does not appear in the other papers. 
For making this comparison we chose the typical value of v=2 corresponding to 7^0.23, 
dc^30Q mi, and 0:^50 mr. Our values of G L in figure 10 are lower than the other estimates, and 
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Figure 10. A comparison of loss in path antenna gain obtained from different scatter theories for a symmetrical 
path for d = 300 miles, M = Sr/6 t =l, s = a //3 = l (antennas on the surface). 



as can be seen in figure 6, they would be even lower if a larger value of v were used or if an 
allowance were made for some degree of anisotropy. The difference between our curve and the 
others is primarily the result of the use of the exponential model atmosphere which is supported 
by ample data [10] and a scattering cross section that is proportional to [2 sin {(a+0)/2}]- 5 . 



The authors acknowledge the assistance of P. L. Kice and H. K. Hansen in the preparation 
of this paper. 
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6. Appendix 1 
The approximation 

1 „ (di-x) a (d a +x) a 



[ 2 sin(^)]W~ d 



~ d b (w 2 +z 2 ) b ' 2 (3-6) 



is usually derived by assuming that a+£ is small [1]. The following method is used to obtain 
this approximation, because it has the advantage that the error is evaluated. Let 

p== R±W. ((M) 



Then, using a well known formula from geometry, 
From the geometry 



7? =V(rf!-x) 2 +2 2 +w 2 (6-3) 



R = ^(d 2 +x) 2 +z 2 +w 2 (6-4) 

and 

«%=(*-,)(*+,) {i+fgp} I/2 { i+ (Sf 2 } 1/2 - (6 - 5) 

(6-5) can be expanded if 2 2 +w 2 <C (di—x) 2 and £ 2 +w 2 <C (rf 2 +#) 2 . This^is generally the case for 
a beyond-t he-horizon propagation path. Then, 

RR = ( dl -x)(d 2 +x)+-2-^ (*-*)(*+*) |+ • • • • (6-6) 

Using (6-3), (6-4), and the first two terms of (6-6), we obtain 

Using (6-2), (6-7), and the first term of (6-6) we obtain 

L. 2sin (^/J (W" ^. {di-xYid^xf 

(RR ) 2 [(R+R Q ) 2 -d 2 Y< 2 ~ dK^+z 2 ) 5 ' 2 ' 

which is the same as (3-6). 
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7. Appendix 2 

The derivation of (3-15) with g mt =l, g mr =l is as follows: 

x-xyidt+xYe- 2 ™ 



«*^»-^{f*£*/>^ 



+ i-a 2 dX )J Z i-J W - (»»+*')'* -\= K ^ K - (Ma) 

where K x and K 2 correspond to the scattering regions to the left and right of the crossover C 
of radio horizon rays in figure 3. 

z lm =Po(d 2 +x). (7-lb) 

z 2m =a (d 1 — x). (7-lc) 

In K 1} let 

3=JW, (7-2) 

and 

So= o^ 1+ ^\ (7 _ 3) 

where 

s=a /p =d 2 ldi. 
In Zs let 

x=-kd 2 (7-4) 

and 

So =^? (1 +S fc). (7-5) 

Then, integrating with respect to w, we obtain 

K °=J^{£ (i -w dk S! z2 ° ex p [? g+w]** 

+sf 1 (l-fc) 3 rfA:r 1 4expr^(l+syi:)lrf3o]-' (7-6) 
where we have set 



Let 



and 



v = 2yh =2ya d x = 2yp d 2 . (4-3) 

£i=iq^ f (1 ~* )3 * f s ° exp [if ft**/ 8 )] rfs °> (7 ~ 7a ) 

L 2 =r4- I \l-k) 3 dk Pag exp [ — (l+«fc)~| eteo- (7-7b) 

i+s Jo Jo L s o J 

We need only consider (7-7a) since (7-7b) can be obtained from (7-7a) by replacing s by 1/s. 
Thus, performing the integration with respect to k, we obtain: 

T l fLi PY „ / ,~-J sz o 3(aa,) 2 . 6(sg ) 3 6(^ ) 4 6 fe ) 4 exp [-^ (ag,)]^ , r7 R . 

L,=— j^exp (- y /z )^ — ^+~ Jz (^ + ~ "^ / d2 °- (7 " 8) 
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Let 

and 
Then 



—=u, 

SO 

_s 4 +7s 3 +2U 2 +35s 



(7-9) 
(7-10) 

(7-11) 



l+ s Jp L u ^ u u 8 J u 8 

f [ ^ ,^+7** S >+7s*+2U* Cl / 6 2 1 1V1 c 



7(1+*') 



where 



+8^(1+1/.)' [i-i+^ 1 -^ + ^-^+^-+4 ^(-„ 2 )]}. 



(7-12) 
(7-13) 



Thus, adding L, and £ s we obtain: 

\8d 

K «= 2lO v W(s+iy ^"1- 720 (^ 7 +l) + 120 (S+7s*+7«+l)v-24 v %,+ iy 

-35^(s+l)]+(s+l)V(6-2i>+^-^)]-(s+l)V[Ei(-p)-Ei(-y[l+l/a]) 
-^i(-j'[l+s])]+se- , ' (1+1 /'' ) [720s 6 -120s 6 (l+s)j»+248 4 (l+s)V-6s 3 (l+s)V 

+2s 2 (l+s)V-.s'(l+.v)V+(l+,s')V] + e-- <1 + s) [720-120(l+.v)»'+24(l+.s') 2 '' 2 

-6(l+s)V4-2(l+s)V-(14-«)V4-(14-s)V]] 



(7-14) 



Figure 1 1 shows 3 ir0 3 (l +s) K /2\Sd plotted versus 77/2 for several values of the path asym- 
metry factor, s. 
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Figuke 11. Dependence of the scattering integral on s and 17. 
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